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ABSTRACT 

Atmospherical mesoscale models can offer unique potentialities to characterize and discriminate potential as- 
tronomical sites. Our team has recently completely validated the Meso-Nh model above Dome C (Lascaux et 
al. 2009, 2010). Using all the measurements of Cjy profiles (15 nights) performed so far at Dome C during 
the winter time (Trinquet et al. 2008) we proved that the model can reconstruct, on rich statistical samples, 
reliable values of all the three most important parameters characterizing the turbulence features of an antarctic 
site: the surface layer thickness, the seeing in the free atmosphere and in the surface layer. Using the same 
Meso-Nh model configuration validated above Dome C, an extended study is now on- going for other sites above 
the antarctic plateau, more precisely South Pole and Dome A. In this contribution we present the most important 
results obtained in the model validation process and the results obtained in the comparison between different 
astronomical sites above the internal plateau. The Meso-Nh model confirms its ability in discriminating between 
different optical turbulence behaviors, and there is evidence that the three sites have different characteristics 
regarding the seeing and the surface layer thickness. We highlight that this study provides the first homogeneous 
estimate, done with comparable statistics, of the optical turbulence developed in the whole 20-22 km above the 
ground at Dome C, South Pole and Dome A. 
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1. INTRODUCTION 

The Internal Antarctic Plateau represents a potential interesting location for astronomical applications. Since 
almost a decade astronomers showed more and more interest towards this region of the Earth thanks to its 
peculiar atmospheric conditions. The extreme cold temperature, the dry atmosphere, the fact that the plateau 
is at more than 2500 m from the sea level, the turbulence seems to develop mainly in a thin surface layer of 
the order of 30-40 m on the top of summits and the seeing above this surface layer assumes values comparable 
to those obtained at mid-latitude sites get this region of the earth very appealing for astronomers. South Pole 
has been the first site equipped with an Observatory in the Internal Antarctic Plateau in which measurements 
of the optical turbulence have been done (Marks et al. 1996, Marks et al. 1999). Fifteen balloons have been 
launched in the winter period and it has been observed that the seeing above a surface layer of ~ 220 m was 
very good (0.37 arcsec). Measurements of the optical turbulence at Dome C are more recent. After the first 
observations done in 2004 with a MASS (Lawrence et al. 2004) , appeared a series of studies done with different 
instrumentation that supported the assessment of the integrated seeing (Aristidi et al. 2005, Aristidi et al. 2009) 
and the vertical distribution of the optical turbulence (Trinquet et al. 2008). 

In this paper we treat a different approach to the site assessment. In this context we are interested in inves- 
tigating the abilities of a mesoscale model (Meso-Nh) in reconstructing correct optical turbulence features above 
different sites of the Internal Antarctic Plateau and its abilities in discriminating the optical turbulence prop- 
erties of different sites. Meso-Nh (Lafore et al., 1998) is a non- hydrostatic mesoscale research model developed 
jointly by the Centre National des Recherches Meteorologiques (CNRM) and the Laboratoire d'Aereologie de 
Toulouse, France. The Astro-Meso-NH package (Masciadri et al. 1999a) has been proved for the first time to be 
able to reconstruct realistic C% profiles above astronomical sites by Masciadri et al. (1999b) and Masciadri et 
al. (2001) and statistically validated later on (Masciadri & Jabouille, 2001, Masciadri et al. 2004, Masciadri et 
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Table 1. Results obtained in Lascaux et al. 2010 that proved the Meso-Nh model reliability above Dome C in reconstructing 
the optical turbulence spatial distribution. Three parameters are estimated: the mean surface turbulent layer (h B ;), the 
seeing in the free atmosphere (sfa) obtained integrating the C% from h s ; up to the end of the atmosphere, the total 
seeing (etot) obtained integrating the C% from the ground up the top of the atmosphere. Beside each parameter is 
reported the associated standard deviation (er) and the statistical error (cr/\/N). 
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Table 2. Geographic coordinates of Dome A, Dome C and South Pole. The altitude is in meter. 



SITE 


LATITUDE 


LONGITUDE 


MESO-NH 


MEASURED 








ALTITUDE (m) 


ALTITUDE (m) 


Dome A* 


80°22'00"S 


077°21'11"E 


4089 


4093 


Dome C** 


75°06'04"S 


123°20'48"E 


3230 


3233 


South Pole 


90°00'00"S 


000°00'00"E 


2746 


2835 



GPS measurement by Dr. X. Cui (private communication). 
GPS measurement by Prof. J. Storey (private communication). 



al. 2006). More recently, it has been statistically validated above Dome C by Lascaux et al. (2009, 2010). The 
most important results obtained in these last papers are summarized in Table [TJ 

Briefly, the observations at Dome C, for a set of 15 winter nights (all the available nights for which is known 
the optical turbulence vertical distribution), gave a mean surface layer thickness h s ^ b s = 35.3 ± 5.1 m. The 
simulated surface layer thickness obtained with the Meso-NH model (h s i :Tnn h — 44.2 ± 6.6 m) is well correlated 
to measurements. The statistical error is of the order of 5-6 m but the standard deviation (a) is of the order 
of 20-25 m. This indicates that the statistic fluctuation of this parameter is intrinsically quite important. The 
median simulated free-atmosphere seeing (e mn h,FA = 0.30 ± 0.17 arcsec) as well as the the median total seeing 
(e mn h,TOT — 1-70 ± 0.21 arcsec) are well correlated to observations, respectively e b s ,FA — 0.3 ± 0.2 arcsec and 
£obs,Tor = 1-6 ± 0.2 arcesc. 

In the context of this paper we consider that the Meso-Nh model is calibrated as shown in Lascaux et al. 
(2010) i.e. it produces optical turbulence features in agreement with observations. We therefore apply the 
Meso-Nh model with the same configuration to other two sites of the plateau: South Pole and Dome A (Table 
0). 

Why these sites ? Dome A is an almost uncontaminated site of the plateau. It is the highest summit of the 
plateau and, for this reason, it is expected to be among the best astronomical sites for astronomical applications. 
The high altitude reduces the whole atmospheric path for light coming from space and above the summit the 
katabatic wind speed is reduced to minima values. Dome A has been proved to have the strongest thermal 
stability (Hagelin et al. 2008) in proximity of the ground due to the coldest temperature. Dome A is a Chinese 
base. In the last few years the Chinese astronomers gave a great impulse to the site characterization showing a 
great interest for building astronomical facilities in this site. Optical turbulence measurements during the winter 
time are not yet available. South Pole is interesting in our study because measurements of optical turbulence 
are available and, at the same time, the site is not located on a summit but above a gently slope. From the 
preliminary measurements done in the past we expect a surface turbulent layer that is thicker than the surface 
layer developed above the other two sites (Dome C and Dome A) due to the ground slope and the consequent 
katabatic winds in proximity of the surface. The three sites form therefore a perfect sample for a benchmark 
test on the model behavior and the model abilities. 

In Section [5] the numerical set-up of the model is presented. In Section [3] results of the complete analysis of 
the three major parameters that characterize the optical turbulence features: surface layer thickness, seeing in 
the free atmosphere i.e. calculated above the surface layer and total seeing are reported. Two different criteria 



Table 3. Meso-NH model configuration. In the second column the horizontal resolution AX, in the third column the 
number of grid points and in the fourth column the horizontal surface covered by the model domain. 
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Domain 1 


25 


120x120 


3000x3000 


Domain 2 


5 


80x80 


400x400 


Domain 3 


1 


80x80 


80x80 



to define the surface layer are used with consequent double treatment. Finally, in Section 0] the results of this 
study are summarized. 

2. NUMERICAL SET-UP 

A detailed description of the Meso-NH mesoscale modeP used for this study is presented inP^ We briefly recall 
here the main characteristics of the numerical configuration: 

• The interactive grid-nesting technique^ is used, with three imbricated domains of increased horizontal 
mcsh-sizes (AX=25 km, 5 km and 1 km, Table [3]). Such a method is used to permit us to achieve the 
best resolution on a small surface but keeping the volumetric domain in which the simulation is done in 
thermodynamic equilibrium with the atmospherical circulation that evolves at large spatial scale on larger 
domains. 

• The vertical grid is the same for all the domains reported in Table [3] The first vertical grid point is at 2 m 
above ground level (a.g.l.). A logarithmic stretched grid up to 3500 m above ground level (a.g.l.) (with 12 
points in the first hundred of meters) is employed. Above 3500 m a.g.l., the vertical resolution is constant 
(AH ~ 600 m). The maximum altitude achieved is around 22 km a.g.l.. 

• All simulations are initialized and forced every 6 hours at synoptic times (00:00, 06:00, 12:00, 18:00) UTC 
by analyses from the European Center for Medium-range Weather Forecasts (ecmwfJJ Note that the 
time at which the simulation starts (UTC) differs for Dome A, Dome C and South Pole. This is done so 
to be able to compare optical turbulence profiles simulated in the same temporal interval with respect to 
the local time (LT). For each night, a mean vertical profile of is computed between the time interval 
(20:00 - 00:00) LT as does in Lascaux et al. (2009, 2010). In Table Hare reported, for each site, the time 
at which the simulation starts and the duration AT of the simulation with respect to the local time. 

• An optimized version of the externalized surface scheme ISBA (Interaction Soil Biosphere Atmosphere) for 
antarctic conditions is employed!^ Such a scheme has been used in Lascaux et al. 2010 and it contributed 
to provide a realistic reconstruction of the optical turbulence near the surface (optical turbulence strength 
and turbulence layer thickness). It is indeed obvious that the most critical part of an atmospherical model 
for this kind of simulations is the scheme that controls the air/ground turbulent fluxes budget. Our ability 
in well reconstructing the surface temperature T s is related to the ability in reconstructing the sensible 
heat flux H that is responsible of the buoyancy- driven turbulence in the surface layer. 

• The Astro-Meso-NH package (Masciadri et al. 1999) implemented in the most recent version of Meso-Nh 
has been used to calculate the optical turbulence and derived astroclimatic parameters. 

As shown inP the best choice for the description of the orography is the RAMP (Radarsat Antarctic Mapping 
Project) Digital Elevation Model (DEM) presented inP instead of the GTOPO30 DEM from the U.S. Geological 
Survey used inP For this study, therefore, the RAMP Digital Elevation Model has been used. The orography 
of each area of interest in this study (Dome C, Dome A, South Pole) is displayed on Fig. [T] All the grid-nested 
domains, from low horizontal resolution (larger mesh-size) to high horizontal resolution (smaller mesh-size) are 

*ECMWF: http://www.ecmwf.int/ 



Table 4. Simulation starting time and time interval chosen for C% computations for the 3 different sets of simulations 
(Dome A, Dome C and South Pole). 





Dome A 


Dome C 


South Pole 


Starting time 


06:00 UTC / 11:00 LT 


00:00 UTC / 08:00 LT 


12:00 UTC / 12:00 LT 


Time interval 
for computations 
(20:00 - 00:00 LT) 


15:00 - 19:00 UTC 


12:00 - 16:00 UTC 


20:00 - 00:00 UTC 



displayed. As can be seen in Fig. [T] (c,f,i) the orography around Dome C and Dome A is more detailed than the 
orography in proximity of the South Pole. This is due to the fact that the procedure to obtain a DEM integrates 
data from many different sources (satellite radar altimetry, airborne surveys, GPS surveys, station-based radar 
sounding...). However the resolution of some areas (typically those that can hardly receive information from the 
satellites) remain poorer than others. The region included in the inner circular polar region (and therefore South 
Pole) fits with this condition and this is the reason why the orography is somehow less detailed than the rest of 
the Internal Antarctic Plateau. Nevertheless, this is a region with no peaks or mountains and with just a regular 
and gently slope. We can therefore reasonably expect that the poorer accuracy in the orography has little or 
minor influence on the results of the numerical simulations done with a mesoscale model such as Meso-NH. 
The same set of 15 winter nights used bj3^ to validate the model above Dome C is investigated in this study 
for the three antarctic sites Dome C, Dome A and South Pole. 



3. OPTICAL TURBULENCE ABOVE DOME C, DOME A, SOUTH POLE 

In this section we investigate and compare the values obtained above the three sites (Dome C, Dome A and 
South Pole) of three parameters that characterize the optical turbulence features above the antarctic plateau: 

• surface layer thickness; 

• free atmosphere seeing from the the surface layer thickness (h s ;) up to the top of the atmosphere; 

• total seeing from the ground up to the top of the atmosphere. We note that this corresponds to ~10 km 
above ground level because the baloons explode at this altitude due to the high pressure and the strong 
wind speed. 

3.1 Optical Turbulence surface layer thickness 

To compute the surface layer thickness for each night, the same method employed irPSl ancP is first used. The 
thickness h s i is defined as the vertical slab containing 90 per cent of the optical turbulence developed inside the 
first kilometer above the ground: 

J8m NV ' < 0.90 (1) 



C m c%{h)dh 



where C\ is the refractive index structure parameter. We remind here that the selection of this criterium (that 
we call criterium A) is motivated by the fact that we intend to compare our calculations with measurements 
done by Trinquet et al. (2008). This criterium has been selected by Trinquet et al. (2008) because the typical 
optical turbulence features above the internal antarctic plateau is characterized by a major bump at the surface 
and a consistent decreasing of the optical turbulence strength in the first tens of meters. The selection of the 
percentage is obviously absolutely arbitrary and it assumes a value only when the same criterium is applied in 
other sites. 

Table [5] reports the computed values of the surface layer thickness for each night at the three sites, as well 
as the mean, the standard deviation (a) and the statistical error {o /\f~N) for the 15 nights. For each night, the 
surface layer thickness is computed from a computed Cjj profile averaged between 20 LT and 00 LT (see Table 
[3] for hours in UT) as done in Lascaux et al. 2010. The calculated mean surface layer thickness above South 
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Figure 1. Orography of three different regions of the internal Antarctic Plateau as seen by the Meson-Nh model (polar 
stereographic projection, grid-nesting configuration), (a), (b) and (c) show the three imbricated domains for the Dome C 
simulations, with horizontal resolution of 25 km, 5 km and 1 km, respectively, (d), (e) and (f) show the three imbricated 
domains for the Dome A simulations, with horizontal resolution of 25 km, 5 km and 1 km, respectively, (g), (h) and (i) 
show the three imbricated domains for the South Pole simulations, with horizontal resolution of 25 km, 5 km and 1 km, 
respectively. The dot labeled 'C indicates the Concordia Station. The dot labeled 'A' indicates Dome A. SP stands for 
South Pole. The altitude is expressed in meter (m). 



Table 5. Mean surface layer thicknesses h s i computed for the 3 sites, for 15 different winter nights using the criterion in 
Eq. [T] Units in meter (m). The mean values are also reported with the associated statistical error a/ '\fN . 
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* These values are computed without taking into account the 
night of the 05/09/05 for Dome C and the night of 07/07/05 
for Dome A (sec text for further explanations). 



Pole (h s /=165 m ± 14.3 m) is more than three time larger than that developed above Dome C (h;,./=43.5 m ± 
5.9 m) and Dome A (h s /=44.2 m ± 8.2 m). This difference is well correlated with previous observations done 
above South Pole. More precisely, observations related to 15 balloons launched during the period (20/6/1995 
- 18/8/1995) indicated h si =220 mP^l Measurements in that paper are done in winter but in a different year 
and different nights. It is not surprising therefore that the matching between calculations and measurements is 
not perfect. Unfortunately the precise dates of nights studied in the paper from Marks et al. (1999) are not 
known. It is therefore not possible to provide a more careful estimate. It is however, remarkable the fact that 
the h s i above South Pole is substantially larger than the h s i above Dome C and Dome A. Also we note that the 
typical thickness calculated above South Pole with a statistical sample of three months bjffSlwas h s / = 102 m. The 
authors used however a different definition of turbulent layer thickness. More precisely, they defined h s ; as the 
elevation (starting from the lowest model level) at which the turbulent kinetic energy contains 1 per cent of the 
turbulent kinetic energy of the lowest model layer. A comparison of this result with our calculations and with 
measurements is therefore meaningless. The same conclusion is valid for the estimates of h s i given at Dome C as 
already explained in Lascaux et al. (2009, 2010). In conclusion, looking at Table O individuals values for each 
nights show a /i s j,5p almost always higher than 100 m, with a maximum close to 300 m (2005 July 25), whereas 
h s i. dc and /i s i,da are always below 100 m. Dome C and Dome A have a comparable surface layer thickness. 
For this sample of 15 nights, h s i,oA is 6.3 m smaller than h s ^oc- We note also that the number of nights for 
which h s i is very small (inferior at 30 m) is more important at Dome A (nine instead of six at Dome C). This 
difference is however not really statistically reliable considering the number of the nights in the sample. For a 
more detailed discrimination between the h s i value at Dome C and Dome A we need a larger statistic. This 
analysis is planned for a forthcoming paper. 

Looking at the results obtained night by night we can note some specific features observe in specific cases. 
Two nights (September 5 at Dome C (h s i = 338.4 m) and July 7 at Dome A (h s i = 529.4 m)) present similar 
characteristics: the surface layer thickness h s ; is well larger than the observed one. In these two cases, however, 
as already explained irF^ for the Dome C site, the large value of h s ; does not mean that a thicker and more 



developed turbulence is present near the ground but it simply means that, in the first kilometer from the ground, 
90 per cent of the turbulence develops in the (0, h s ;) range. In both cases the model reconstructs the total 
seeing on the whole 20 km much weaker than what has been observed and more uniformly distributed and, 
consequently, the criterium (EqJTJ) provides us a much larger value of h s ;. In both cases, when we look at the 
vertical distribution of the C^r , we observe that the turbulence is concentrated well below 20 m in a very thin 
surface layer with a very weak total seeing (see next section and table [5]). The case of 5 September at Dome A, is 
however a case in which the model reconstructed a surface turbulent layer thicker than what has been observed. 

In order to compare our calculations and results with those obtained by^we applied also a different criterium 
(criterium B) based on the analysis of the vertical profile of turbulent kinetic energy (TKE) instead of the vertical 
profile of C^f . The TKE is certainly an ingredient from which the optical turbulence depends on and it represents 
the dynamic turbulent energy. However, it is known (Masciadri & Jabouille 2001) that the C]^ depends also on 
the gradient of the potential temperature and moreover, it is absolutely arbitrary the selection of the value of 
percentage of the turbulent kinetic energy (1%, 10%, other...) used as a threshold. This method is not therefore 
useful to quantify the absolute value of h s j to be compared to measurements provided by Trinquet et al. (2008) 
and Marks et al. (1999). It can be, eventually, useful for relative comparisons between different sites or to 
compare our calculations with calculations provided by Swain & Gallee (2006). 

Using this method (Tabled, the surface layer height is determined as the elevation at which the TKE is X% 
of the lowest elevation value. We calculated the h s ; for X = 1 (Table O and X = 10 (Table [7]). X=l is the case 
treated byP^l For each simulation, we first compute the average of the TKE profile for the night between 20 LT 
and 00 LT. While the average of the profile is calculated with a 2 minutes rate sample, the average of the 
TKE is calculated with 5 profiles, available at each hour (20, 21, 22, 23, 24, 00) LT. This gives us an averaged 
vertical profile of TKE characteristic of the considered night. 

The computation of the surface layer thickness is then performed using this averaged TKE profile. It has 
been observed that, when the night presents only low dynamic turbulence (with a very low averaged TKE at 
the lowest elevation level), it is very hard to retrieve a surface layer height using this criterium. This means that 
the turbulence is so weak that we are at the limit of necessary turbulent kinetic energy to resolve the turbulence 
itself. For these nights (indicated with an asterisk in Table it could happen that we calculated the average on 
a number of estimates minor than 5 (as for all the other cases). The results are reported in Table El 

Table [6] shows that results obtained with the criterium of the TKE are similar to those obtained with the 
criterium described in Eq|T] Table [7] provides, however smaller values of h s ; above all the three sites. We treat 
the case (X = 10) to show that, tuning the value of the percentage, it is possible to find different values of h s ;. 
This means that estimates are useful only if they are compared to measurements using the same criteria. 

To conclude, both criteria (A and B with X=l) give similar mean h s i values for all the 3 sites for this limited 
set of nights. Evaluating the surface layer thickness over a more extended set of nights should be the next step. 
It would permit us to compute more reliable and robust statistical estimates for h s i over the 3 antarctic sites. 
This also means that our estimate (h s ; = 165 m) above South Pole is better correlated to measurements (h s ;=220 
m) that obtained bjM (h s j = 102 m) at the same site. 

3.2 Seeing in the free atmosphere and seeing in the whole atmosphere 

Table \E\ shows the simulated total seeing (etot) and free-atmosphere seeing {efa) for each night and each sites 
(Dome C, Dome A and South Pole). We define the free atmosphere as the portion of the atmosphere extended 
from the mean h s i reported in Table [5] up to 10 km. At this height the balloon explode in general during the 
winter and we need to limit the integration of the at such a low value to be able to compare our results 
with those obtained with measurements. The median values of the seeing as well as the standard deviation (a) 
and the statistical error (cr/yN) are reported. As expected the total seeing is a little bit stronger at Dome A 
{etot,da = 2.37 ± 0.27 arcsec) than at Dome C (stot,dc — 1-70 ± 0.21 arcsec) or South Pole (etot,sp = 1-82 
± 0.23). The total seeing is very well correlated with measurements at Dome C (Lascaux et al. 2010 - £TOT : obs = 
1.6 arcsec) and at South Pole (Marks et al. 1999 - £TOT,obs = 1-86 arsec) getting the estimate at Dome A highly 
reliable. The minimum median free-atmosphere seeing at Dome A is Sfa.da — 0.23 ± 0.28 arcsec, at Dome C, 
£fa,dc = 0.30 ± 0.17 arcsec and at South Pole, Sfa.sp = 0.36 ± 0.11 arcsec. The seeing in the free atmosphere 



Table 6. Mean surface layer thicknesses h s i computed for the 3 sites, for the same set of nights shown in Table [S] but 
computed with a different criterion. The surface layer height is determined as the elevation at which the averaged TKE 
between 20 LT and 00 LT for each nigh is 1% of the averaged lowest elevation value. Units in meter (m). The mean 
values are also reported with the associated statistical error a/y/N. 
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* These values are computed using a number 
of profiles minor than five. 



Table 7. Mean surface layer thicknesses h s i computed for the 3 sites, for the same set of nights shown in Table [5] but 
computed with a different criterion. The surface layer height is determined as the elevation at which the averaged TKE 
between 20 LT and 00 LT for each nigh is 10% of the averaged lowest elevation value. Units in meter (m). The mean 
values are also reported with the associated statistical error g/\/N. 
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Table 8. Total seeing ETOT=£[8m,,h t o ] anc ^ seeing in the free atmosphere £FA~£[h 3l ,h to ] calculated for the 15 nights and 
averaged in the temporal range 20-00 LT . See the text for the definition of h s ; and ht op . In the second column are 
reported the observed values, in the third and fourth columns the simulated values obtained with high and low horizontal 
resolution respectively. Units in arcsec. 





DOME A 


DOME C 


SOUTH POLE 


Date 


£fa/£tot 


£fa/£tot 


£fa/ £tot 




(fc s( =37.9ra) 


(h sl =44. 2m) 


(h gl =165m) 



04/07/05 


2.55 


/ 


3 


37 





22 


/ 


2.28 


0.40 


/ 


1 


67 


07/07/05 


0.20 


/ 





24 





28 


/ 


1.91 


0.31 


/ 





70 


11/07/05 


0.23 


/ 


2 


78 


1 


61 


/ 


1.81 


0.47 


/ 


1 


96 


18/07/05 


0.21 


/ 


2 


73 





80 


/ 


1.94 


1.46 


/ 


2 


28 


21/07/05 


0.21 


/ 


1 


95 





86 


/ 


1.27 


0.31 


/ 


1 


71 


25/07/05 


0.22 


/ 


1 


55 





25 


/ 


0.85 


0.32 


/ 





76 


01/08/05 


0.22 


/ 


1 


78 





22 


/ 


2.27 


0.52 


/ 


1 


78 


08/08/05 


1.45 


/ 


2 


42 





35 


/ 


1.70 


0.28 


/ 


1 


69 


12/08/05 


0.23 


/ 


2 


37 





23 


/ 


0.99 


0.29 


/ 


1 


82 


29/08/05 


0.23 


/ 


1 


83 


2 


29 


/ 


2.47 


1.55 


/ 


2 


11 


02/09/05 


0.22 


/ 


1 


7G 


1 


16 


/ 


1.54 


0.81 


/ 


3 


56 


05/09/05 


3.21 


/ 


3 


3G 





30 


/ 


0.52 


0.31 


/ 


2 


98 


07/09/05 


2.43 


/ 


3 


49 


1 


69 


/ 


3.73 


0.31 


/ 


1 


41 


16/09/05 


1.11 


/ 


4 


GO 





21 


/ 


1.57 


0.99 


/ 


3 


96 


21/09/05 


0.20 


/ 


2 


30 





26 


/ 


1.63 


0.36 


/ 


2 


32 



Median 0.23 / 2.37 0.30 / 1.70 0.36 / 1.82 
a 1.08 / 1.03 0.67 / 0.77 0.43 / 0.90 
a/VN 0.28 / 0.27 0.17 / 0.21 0.11 / 0.23 



H (km) H (m) H (m) 




Figure 2. Median C% profiles simulated with the Meso-NH mesoscale model at Dome C (black), Dome A (blue) and 
South Pole (red). Left: from the ground up to 20 km. Middle: from the ground up to 1 km. Right: from the ground up 
to 200 m. Units are m" 2/3 . 



is very well correlated with measurements at Dome C (Lascaux et al. 2010 - £FA,obs = 0.30 arcsec) and at South 
Pole (Marks et al. 1999 - EFA,obs = 0.37 arsec) getting again very reliable the method (Meso-NH model) as well 
as the estimates at Dome A. What is remarkable is that, even if h s i t DA < ^s/.dc < h s i,sp> Dome A is the site 
with the lowest free-atmosphere seeing Efa- This means that at Dome A as well as at Dome C the turbulence 
is concentrated inside the first tenths of meters from the ground. Moreover, the turbulence in the surface layer 
is stronger at Dome A than at Dome C. This can be explained with the stronger thermal stability of Dome A 
near the ground. Our results match, therefore, with predictions we did in Hagelin et al. (2008) studying only 
features of the meteorological parameters. At South Pole, however, the C]^ vertical distribution decreases in 
a less abrupt way. and it is spread over hundredth of meters from the ground, instead of tenths of meters like 
for Dome A or Dome C. As a result the total seeing is also weaker than above Dome C and Dome A. Such a 
behavior is evidenced in Figure [U which displays the median vertical C^, profiles over the 3 sites. 

4. CONCLUSION 

In this study the mesoscale model Meso-NH was used to perform forecasts of optical turbulence (evolutions of 
Cff profiles) for 15 winter nights at three different antarctic sites: Dome A, Dome C and South Pole. The model 
has been used with the same configuration previously validated at Dome C (Lascaux et al. 2010) and simulations 
of the same 15 nights have been performed above the three sites. The idea behind our approach is that once 
validated above Dome C, the model can be used above two other site above the internal antarctic plateau to 
discriminate optical turbulence features typical of other sites. This should show the potentiality of the numerical 
tool in the context of the site selection and characterization in astronomy. South Pole has been chosen because in 
the past some measurements of the optical turbulence have been done and this can represent a useful constraint 
for the model itself. For Dome A there are not at present time measurements of the optical turbulence and 
this study provides therefore the first estimates ever done of the optical turbulence above this site. We test this 
approach above the antarctic plateau because this region is particularly simple from the morphologic point of 
view and certainly simpler than typical mid-latitude astronomical sites. No major mountain chains are present 
and the local surface circulations is mainly addressed by the energy budget air/ground transfer, the polar vortex 
circulation at synoptic scale and the katabatic winds generated by gravity effects on gently slopes due to the 
cold temperature of the iced surface. 

The main results we obtained are: 

• The Meso-Nh model achieves to reconstruct the three most important parameters used to characterize 
the optical turbulence: the turbulent surface layer thickness, the seeing in the free atmosphere and in the 
surface layer for the three selected sites: Dome C, Dome A and South Pole showing results in agreement 
with expectations. Measurements taken at Dome C and South Pole corresponds to balloons launched 
during 15 nights, in both cases. The statistic is not very large but reliable for a first significant result. The 
selected nights correspond to the 15 nights for which measurements of the Dome C are available. 

• Dome C and Dome A present a very thin surface layer size (Ii^da — 37.9 ± 8.1 m and h s i t DC = 44.2 
± 6.6 m) while South Pole surface layer is much thicker (/i s ;,sp = 165 ± 17.4 m). All the estimates are 
well correlated with measurements. To better discriminate between Dome A and Dome C surface layer 
thickness a richer statistic is necessary. An on-going study has started addressing this issue. 

• Dome A is the site with the strongest total seeing (2.37 ± 0.27 arcsec) with respect to Dome C (etot,dc 
= 1.70 ± 0.21 arcsec) and South Pole (etot.sp — 1-82 ± 0.23 arcsec). This is explained by the stronger 
thermal stability near the ground with respect to other two sites that cause large values of the optical 
turbulence in the thin surface layer. 

• All the three sites show a very weak seeing in the free atmosphere i.e. above the correspondent mean h s ;: 
efa.da = 0.23 ± 0.28 arcsec at Dome A, Efa.dc = 0.30 ± 0.17 arcsec at Dome C and Efa.sp — 0.36 ± 
0.11 arcsec at South Pole. Dome A show the weakest seeing in the free atmosphere. 

• Both, the total seeing and the seeing in the free atmosphere calculated by Meso-Nh, are very well correlated 
with measurements at Dome C and South Pole getting the predictions done at Dome A highly reliable. 



• Dealing on the criteria used to define the surface layer thickness, we proved that, at least on the sample 
of 15 nights investigated, the criterium defined by EqfT] (criterium A) and the criterium using the vertical 
profile of the turbulent kinetic energy (TKE) taking h s ; as the height at which the value of the TKE is 
less than 1% of the TKE at the lowest level near the ground (criterium B) provide more or less the same 
results. 

• The mean h s ; we estimate at Dome C (h s ;=44.2) is slightly thicker than what found bjff^ (h s ;=27.7 m) 
with comparable discrepancy from measurements (h s i = 35.3 ± 5.1 m). The h s ; we estimate at South 
Pole (h s isp = 165 ± 17.4 m) is thicker than what estimated bjff^ (/i s ;,sp = 102) but better correlated to 
measurements (h s i t sp = 220 m) than what found byP^I The h s / we estimate at Dome A {h s i^A — 37.9 
± 8.1 m) is somehow thicker than what estimated by^l (Ji s i_£ia = 18 m). It is however important to note 
that the standard deviation of h s i is of the order of h s ; itself or even larger. The statistic error a/ yf{N) 
is of the order of ~ 10 m. We think therefore that at present there are no major differences in our results 
with respect tcE^with expection of the fact that we proved that, with our model, the horizontal resolution 
of 1 km provides better results than a resolution of 100 km (G^). 

All these results deserve now a confirmation provided by an analysis done with a richer statistical sample. 
Besides, we can state that all major expectations concerning the typical features of the optical turbulence above 
South Pole, Dome C and Dome A have been confirmed by this study. We highlight that this paper provides the 
first estimate of the optical turbulence extended on the whole 20 km above the Internal Antarctic Plateau. 
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